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Narrow Line Seyfert 1 galaxies (NLSls) are generally considered peculiar objects 
among the broad class of Type 1 active galactic nuclei, due to the relatively small 
width of the broad lines, strong X-ray variability, soft X-ray continua, weak [O in] , and 
strong Fen line intensities. The mass Mbh of the central massive black hole (MBH) is 
claimed to be lighter than expected from known MBH-host galaxy scaling relations, 
while the accretion rate onto the MBH larger than the average value appropriate to 
Seyfert 1 galaxies. In this Letter, we show that NLS1 peculiar Mbh an d L/LEdd turn 
out to be fairly standard, provided that the broad line region is allowed to have a 
disc-like, rather than isotropic, geometry. Assuming that NLSls are rather "normal" 
Seyfert 1 objects seen along the disc axis, we could estimate the typical inclination 
angles from the fraction of Seyfert 1 classified as NLSls, and compute the geometrical 
factor relating the observed FWHM of broad lines to the virial mass of the MBH. We 
show that the geometrical factor can fully account for the "black hole mass deficit" 
observed in NLSls, and that L/L^dd is (on average) comparable to the value of the 
more common broad line Seyfert 1 galaxies. 
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1 INTRODUCTION 

Seyfert 1 galaxies (Syls) are often divided into two distinct 
classes, namely Broad Line Syls (BLSls), whose H/3 line has 
FWHM > 2000 km/s (hence, as standard Type 1 AGN), 
and Narrow Line Syls (NLSls), with lower velocities (e.g., 
Goodrich, 1989). NLSls are a minority, ~ 15% of all the 
Syls, according to the optical spectroscopic classification of 
the SDSS general field (Williams Pogge & Mathur, 2002), 
the fraction depending on the AGN luminosity (with a peak 
at M g i ~ —22), and on the radio loudness (radio loud NLSls 
account only for ~ 7% of the class, Komossa et al., 2006, 
but it is still debated if the NLSls can be considered a pe- 
culiar radio-quiet sub-class among Syls, see e.g. Sulentic et 
al. 2007; Doi et al. 2007). NLSls also show weak [Om] and 
strong Fell emission line (Osterbrock & Pogge, 1985), strong 
variability, and a softer than usual X-ray continuum (Boiler 
Brandt & Fink, 1996; Grupe et al., 1999). 

Grupe & Mathur (2004) found that NLSls have, on av- 
erage, lower Mbh than expected from AfeH-host galaxy re- 
lations such as Mbh-ct* (see Tremaine et al., 2002, and refer- 
ences therein), while BLS1 A/bh are in fairly good agreement 
to the same relation. The estimated low values of Mbh lead 
to an average Eddington ratio L/LEdd for the NLS1 pop- 



ulation which is almost an order of magnitude larger than 
the average value of BLSls (L/LEdd — 1 to be compared to 
~ 0.1, Grupe, 2004). Further evidence of low Mbh in NLSls 
comes from the observed rapid X-ray variability (see., e.g., 
Green, McHardy & Lehto 1993, and Hayashida 2000). 

Such results were interpreted as indications of a peculiar 
role of NLSls within the framework of the cosmic evolution 
of MBHs and of their hosts. In a MBH-galaxy co-evolution 
scenario, NLSls are thought to be still on their way to 
reach the A/bh-o"* relation, i.e., their (comparatively) small 
MBHs are highly accreting in already formed bulges. Re- 
cently Botte et al. (2005) and Komossa and Xu (2007) came 
to the conclusion that NLSls have indeed smaller masses 
and higher L/L^dd than BLS1, nevertheless they both do 
follow the M — a, relation for quiescent galaxies. The au- 
thors argued that the customarily used [O ill] line is not a 
reliable surrogate for the stellar velocity dispersion a*. 

The Grupe and Mathur's results and interpretation 
have been recently confirmed and supported by several 
other groups, see, e.g., Zhou et al. (2006) and Ryan et al. 
(2007). Ryan et al. (2007) pointed out that IR-based mass 
measurements might be unreliable because of the extra IR 
contribute from the circum-nuclear star-forming regions in 
NLSls. Notwithstanding, they suggested that this contami- 
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nation can not significantly affect their data, and thus is in- 
sufficient to account for the MBH mass deficit. In the afore- 
mentioned papers, Mbh was computed as 



Mbh = 



G 



(1) 



where -Rblr is the broad line region (BLR) scale radius, and 
WBLR the typical velocity of BLR clouds. 7?blr is found by 
means of the reverberation mapping technique (Blandford & 
McKee, 1982), or by exploiting statistical i?BLR-luminosity 
relations (see Kaspi et al., 2000, 2005 and 2007); hblr can 
be inferred from the H/3 width as 



«blr = / ■ FWHM, 



(2) 



where the FWHM refers only to the broad component of 
the line, and / is a fudge factor which depends upon the 
assumed BLR model. For an isotropic velocity distribution, 
as generally assumed, / = v3/2. 

Labita et al. (2006) and Decarli et al. (in preparation) 
found that in QSOs an isotropic BLR fails to reproduce the 
observed line widths and shapes, and a disc model should be 
preferred. A disc-like geometry for the BLR has been pro- 
posed by several authors in the past (e.g., Wills & Browne 
1986; Vestergaard, Wilkes & Barthel 2000; Bian & Zhao, 
2004). In this picture, the observed small FWHM of NLS1 
broad lines are ascribed to a small viewing angle with respect 
to the disc axis, and no evolutionary difference is invoked 
whatsoever. 

In this Letter; we adopt the disc-like model for the BLR 
of Seyfert galaxies. We use the observed frequency of NLSls 
to estimate their typical viewing angle, and then compute 
the appropriate geometrical factor /. Using eq. [T] we will 
show that the new estimates of Mbh for NLSls nicely agree 
with the standard Mbh-ct. relation. In turn, the accretion 
rate of the class is found to be similar to that of BLSls. 



2 MODEL AND RESULTS 

We model the BLR as a thin disc, and define i? as the angle 
between the line of sight and the normal to the disc plane. 
The FWHM is a measure of the velocity projected along the 
line of sight. In the assumption of a 2-D, keplerian BLR, the 
observed FWHM is correlated to the rotational velocity of 
the disc as following: 



FWHM = 2u kcp sin??, 



(3) 



where v^ ep is the keplerian velocity of the disc-like BLR. In 
this model the differences between the FWHM of NLSls and 
BLSls depend only on #, so that the Syls observed nearly 
face-on are classified as NLSls, while the ones observed at 
higher angles are considered BLSls. As mentioned in the 
introduction, the fraction of NLSls we consider is ~ 15%. 
In our unification scheme, the relative fraction -Rnlsi is re- 
lated to the maximum inclination angle of the subclass i? C r 

as i?NLSl = (l-COS$ C r)/(l-COS#max), where 1?max ~ 40° is 

the maximum inclination angle for Type-I AGNs in the uni- 
fication model (e.g., Antonucci & Miller, 1985; Antonucci, 
1993; Storchi-Bergmann, Mulchaey & Wilson, 1993). 

Some authors suggested that the BLR can not be com- 
pletely flat (see, e.g., Collin et al., 2006). Alternatively, discs 
may have a finite half thickness (H), or a "flared" profile 
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Figure 1. The dependence of i? cr (red, solid line) and /nlsi on 
i?max- The blue, dashed line and the magenta, dotted line refer to 
values /nlsi calculated assuming H/R = 0.1 and respectively. 



(with H increasing more than linearly with the disc radius 
R, see, e.g., Dumont & Collin-Souffrin, 1990). Other mod- 
els proposed include warped discs (Wijers & Pringle, 1999), 
and the superposition of discs and wind components (Mur- 
ray & Chiang, 1995, 1998; Elvis, 2000; Proga & Kallman, 
2004). In this Letter we employ the simplest model, i.e., a 
disc with finite thickness, a choice minimizing the number 
of required parameters. As it will be shown in the following, 
this minimal set-up can resolve the apparent dichotomy be- 
tween NLSls and BLSls. 

In a finite thickness disc model for the BLR, the geo- 
metrical factor /, as defined in equation [2] is related to the 
inclination angle $ of the disc as 



(|) 2 + sin^ 



(4) 



The ratio H/R is related to the relative importance of 
isotropic (e.g. turbulent) vs rotational motions. The aver- 
age geometrical factor for each class, /nlsi and /blsi, is 
computed by averaging equation [4] over the relevant solid 
angle (0 < d < # cr for NLSls, # cr < ■& < # max for BLSls). 

Fig. [T]shows the dependence of i9 cr and /nlsi on i? max , 
with 35° < i? ma jt ^ 50°. The critical angle ranges between 
13° and 19°, while /nlsi is found between ~ 3 and 4.5 in 
the limit H/R = 0, and between ~ 2.2 and 2.9 for H/R = 
0.1. We also find 0.9 < /blsi < 1.2 independently of < 
(H/R) < 0.1. 

We adopt a fiducial value i) ma , = 40°, leading to 
/nlsi ^ 3.8 and ~ 2.6 for H/R = and H/R = 0.1, re- 
spectively, and /blsi — 1- 

The new estimates of the geometrical factor allow us 
to reconsider the values of Mbh for the sample of Syls pre- 
sented in Grupe & Mathur (2004), who instead employed 
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Figure 2. Upper panel — the cumulative fraction distribution of 
Grupe & Mathur (2004) sample as a function of Mbh- The black, 
solid line refers to BLSls, after applying our correction. The red, 
dot-dashed line refers to NLSls with /nlsi = V^/2. The blue, 
dashed and the magenta, dotted lines refer to NLSls assuming 
H/R = 0.1 and 0, respectively. Lower panel - the distribution 
of Grupe & Mathur (2004) sample on the Mbh - c* plane. Black, 
empty circles refer to BLSls, when /blsi — 1 is adopted. Red, 
filled squares are NLSI values, using /nlsi = v / 3/2. Blue filled 
triangles refer to the NLSls after the correction described in the 
text, assuming H/R = 0.1. The arrows highlight that the values 
of ct„ for NLSls have to be considered upper limits, as discussed 
in the text. The Tremaine et al. (2002) relation is also plotted for 
comparison. 



a fixed / = \/3/2 for all objects. Our results are shown 
in Fig. [2] In the upper panel the blue long-dashed (ma- 
genta dotted) line refers to the corrected A/bh of NLSls 
for H/R = 0.1 (H/R = 0). NLSI black hole masses are 
increased by ~ 0.84 (~ 1.16) dex, while BLSI black hole 
masses by a mere ~ 0.05 (~ 0.07) dex, with respect to the 
Grupe & Mathur values. The mass distributions for the two 
classes are now remarkably similar, without any significant 
difference between NLSls and BLSls. 

The lower panel of Fig. [2] shows the BLSI and NLSI 
populations in the Mbh-o* plane. The black, empty circles 
refer to BLSls, assuming /blsi — 1- The red, solid squares 
are NLSls for /nlsi = V3/2, while the blue solid triangles 
refer to the NLSls after the correction described in the text 
is applied, assuming H/R — 0.1. The estimates of <r* are 
from Grupe & Mathur (2004), and are derived from [Oiii] 
line width. It should be noted that, as Botte et al. (2005) 
and Komossa & Xu (2007) pointed out, the [O ill] surrogate 
poorly correlates with er* measured from stellar absorption 
lines, so that the plotted a, values have to be considered 
upper limits, as indicated by the arrows. This caveat is par- 
ticularly important for X-ray selected samples, as the one 
used here (Marziani et al., 2003), as wind components to 
[O ill] lines may be significant. 
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Figure 3. Cumulative fractions of NLSls and BLSls as a function 
of the Eddington ratio. The color/line-type code is the same as 
in Fig. 2, upper panel. 



We can now estimate the corrected Eddington ratio for 
the same sample (Grupe, 2004). The cumulative fractions 
of NLSls and BLSls vs L/L^dd are shown in Fig. [3] Not 
surprisingly, having comparable luminosities, and now, com- 
parable masses, NLSls and BLSls radiate at the same Ed- 
dington ratio. This result supports the pole-on orientation 
model for NLSls. 



3 DISCUSSION AND CONCLUSIONS 

In this Letter, we have assessed the claimed peculiarity of 
Narrow Line Seyfert 1 galaxies within the framework of cos- 
mic evolution of massive black holes, and their host bulges. 
Indeed, the optical properties of NLSls, their X-ray fast 
variability and the faintness of their bulges can be accounted 
for if one admits lower black hole masses and higher accre- 
tion rates (in Eddington units) than standard Broad Line 
Seyfert 1 galaxies (BLSls), placing NLSls in an early evolu- 
tionary stage (Grupe & Mathur, 2004; Grupe, 2004; Botte 
et al., 2004; Zhou et al., 2006; Ryan et al., 2007). If this 
is true, by observing local NLSls we can have hints of the 
infancy of the ubiquitous population of super-massive black 
holes. 

We have explored an alternative explanation to the nar- 
rowness of H/3 lines in NLSls, namely, pole-on orientation 
of a disc-like broad line region. If BLSls and NLSls differ 
only by the observation angle of the BLR disc, the frequency 
of NLSls among the Syl class gives the limiting viewing 
angle of NLSls. Then, assuming H/R < 0.1 for the disc, 
we computed corrected geometrical factors linking the ob- 
served FWHM to Mbh, and found /nlsi ^ 2 and /blsi — 1, 
in agreement with recent estimates given by Labita et al. 
(2006). The idea of a disc-like BLR is not new (e.g., Wills 
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& Browne 1986; Vestergaard, Wilkes & Barthel 2000; Bian & 
Zhao, 2004), but for the first time, by re-calculating masses 
and Eddington ratios for a sample of Syls, we found that 
mass and luminosity functions are similar for NLSls and 
BLSls. In a sense, we can say that all Syls are normal, but 
some are more "normal" than others. 

We note that, though NLSls seem to lie in the same 
region of the Mbh-o* plane, the adopted ct* values can be 
largely over-estimated (Komossa & Xu, 2007; Mullaney & 
Ward, 2007), and then firm conclusions on the Mbh-o* issue 
can not be drawn at this stage. 

Can a simple orientation model, as the one we adopted 
here, explain the unique observed properties of NLSls? 
NLSls differ from standard Syls not just in the width of 
optical lines, but, more noticeably, in what are their X-ray 
properties, both spectral and temporal. The X-ray emission 
of NLSls has been studied and discussed in great details by, 
among others, Boiler et al. (1996), using ROSAT data, and 
by Brandt, Mathur & Elvis (1997) using ASCA data. NLSls 
have generally both soft and hard X-ray spectra which are 
steeper than normal Syls, and show large amplitude, rapid 
variability. Boiler et al. (1996) showed how different models, 
invoking extreme values of one or more of the followings: 
pole-on orientation, black hole mass, accretion rate, warm 
absorption, BLR thickness, all explain some aspects of the 
complex NLS1 soft X-ray phenomenology, but, still, all ap- 
pear to have drawbacks. 

If pole-on orientation has to be the main cause of the 
uniqueness of the X-ray features of NLSls, then a neces- 
sary condition is that the hard power-law emission is not 
intrinsically isotropic, e.g., a thermal extended corona (as 
in Haardt & Maraschi, 1991; 1993) is not a viable option. 
Models in which the X-rays of type I radio quiet AGNs are 
funneled or beamed have been proposed by several authors 
(e.g., Madau, 1988; Henri & Petrucci, 1997; Malzac et al., 
1998; Ghisellini, Haardt & Matt, 2004). For example, Ghis- 
ellini et al. (2004) showed that an aborted jet model, in 
which X-rays are produced by dissipation of kinetic energy 
of colliding blobs launched along the MBH rotation axis, 
can explain the steep and highly variable X-ray power law. 
The model, in its existing formulation, does not allow clear 
predictions of spectral and temporal features other than in 
the X-rays. To assess its relevance for NLSls would require 
a much more detailed modeling. In particular, the peculiar 
Fell and [Oni] properties must be accounted for. 

The statistics of radio-loud NLSls is low. In several 
works the existence of differences in the radio properties 
between NLSls and BLSls has been discussed (see, e.g. Ko- 
mossa et al. 2006; Zhou et al. 2006; Sulentic et al. 2007; 
Doi et al. 2007). Doi et al. (2007) suggested that ~ 50 % of 
radio- loud NLSls are likely associated with jets with high 
brightness temperatures, requiring Doppler boosting. This 
interpretation supports the pole-on orientation model (for 
a different point of view see Komossa et al. 2006). 

Our results, if confirmed, indicate that a population of 
accreting, undermassive MBHs (with respect to the Mbh~ 
cr* relation) has to be found yet. This may suggest that 
the Mbh-c* relation was established long ago, during the 
MBH accretion episodes following the first major mergers 
of the host galaxies. Moreover, Komossa & Xu (2007) found 
that NLSls do follow the Mbh-o"* relation of non-active 
galaxies, but still they have smaller Mbh and larger L/Z/Edd 



than BLSls. If this is the case, then a* of the host bulges 
of NLS1 needs to evolve accordingly in order to preserve 
the Mbh-c* relation, or, alternatively, NLSls are the low 
mass extension of BLSls, and the NLS1 high L/I/Edd is a 
short-lived phenomenon. We note here that the interpreta- 
tion of Komossa & Xu (2007), as well as the one of Grupe 
& Mathur (2004), implies that Mbh and L/L^aa, in prin- 
ciple independent quantities, somehow conspire to produce 
comparable luminosities as observed in NLSls and BLSls. 
Applying our correction to the Mbh as well as the one to 
the <r* proposed by Komossa & Xu (2007), the NLSls would 
be even off-setted towards higher masses with respect to the 
Mbh-c* relation. 

There are however two possible problems with the 
pole-on orientation model. First, according to the orienta- 
tion model, the polarization properties of broad emission 
lines should depend on the inclination angle, in the sense 
that nearly pole-on Syls should not be polarized. How- 
ever, Smith et al. (2004) found polarized broad lines in few 
NLSls, and traces of broad Ha polarization were also found 
by Goodrich (1989) in 6 out of 17 NLSls. A second issue is 
discussed by Punsly (2007), who finds larger line broaden- 
ing in face-on quasars, possibly due to large isotropic gas 
velocities or winds. 

In conclusion, we found that orientation effects can ac- 
count for the different optical properties of NLSls compared 
to the more common BLSls. The model is particularly ap- 
pealing, as it naturally sets masses and accretion rates of 
NLS1 to fairly standard values. To validate this interpreta- 
tion, orientation must be able to explain the extreme X-ray 
properties of NLS1. Jetted models for radio quiet AGNs may 
be promising in this, and we urge a detailed, critical com- 
parison of such models with the bulk of NLS1 data. 
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